Maximally polarized states for quantum light fields 
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degree of polarization of a quantum state can be defined as its Hilbert-Schmidt distance to the 
set of unpolarized states. We demonstrate that the states optimizing this degree for a fixed average 
number of photons N present a fairly symmetric, parabolic photon statistics, with a variance scaling 
a.s N 2 . Although no standard optical process yields such a statistics, we show that, to an excellent 
approximation, a highly squeezed vacuum can be considered as maximally polarized. 
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A number of key concepts in quantum optics can be 
concisely quantified in terms of distance measures. The 
notions of nonclassicality 0, entanglement Q, informa- 
tion , localization Q , or polarization [f| , to cite only a 
few relevant examples, have been systematically formu- 
lated within this framework. The rationale behind this 
is quite clear: once we have identified a convex set with 
the desired physical properties (classicality, separability, 
etc.), the distance determines the distinguishability of a 
state with respect to that set 0. Apart from its con- 
ceptual simplicity, this procedure avoids many undesired 
problems that can arise in more standard approaches. 

Irrespective of our particular choice for the distance, 
a natural question emerges: what states maximize the 
corresponding measure. A good deal of effort has been 
devoted to characterize maximally nonclassical or entan- 
gled states. However, as far as we know, maximally po- 
larized states have been not considered thus far, except 
for some trivial cases. It is precisely the purpose of this 
Letter to fill this gap, providing a complete description 
of such states, as well as feasible experimental schemes 
for their generation. 

Let us start by briefly recalling some basic concepts 
about quantum polarization. We assume a monochro- 
matic plane wave propagating in the z direction, whose 
electric field lies in the xy plane. Under these condi- 
tions, the field can be fully represented by two complex 
amplitude operators, denoted by an and ay when using 
the basis of linear (horizontal and vertical) polarizations. 
They obey the standard bosonic commutation relations 
[aj,ajy = Sjk, with j,k £ {H, V}. The Stokes operators 
are then introduced as the quantum counterparts of the 
classical variables 7] , namely 



So = a' H aH + a v dv 



Si = i(a H ay - a* H a v ) , 



Si = a H av + a-v a H 



(1) 



S 3 = a} H a H - OyCiv , 



and their mean values are precisely the Stokes parameters 
((So),(S)), where S T = (Si, Si, S3) and the superscript 



T indicates the transpose. One immediately gets that the 
Stokes operators satisfy the commutation relations of the 
algebra su(2): [Si, 52] = 2iSa, and cyclic permutations. 
Moreover, since [S, So] = 0, each energy manifold (with 
a fixed number of photons) can be treated separately. To 
bring out this point more clearly, it is advantageous to 
relabel the standard two-mode Fock basis \n)jj <S> \m)v 
in the form 



\N, k) = \k) H ® \N - k) v , 



k = 0,l, 



,N. (2) 



These states span an invariant subspace of dimension 
N + 1, and the operators S act therein as an angular 
momentum N/2. 

Any (linear) polarization transformation is generated 
by the Stokes operators Q. However, So induces only a 
common phase shift that does not change the polariza- 
tion state and can thus be omitted. Therefore, we restrict 
ourselves to the SU(2) transformations, generated by S. 
Since Si is related to S2 and S3 by the commutation rela- 
tions, only these two generators suffice. It is well known 
that S2 generates rotations around the direction of prop- 
agation, whereas S3 represents differential phase shifts 
between the modes. It follows then that any polarization 
transformation can be realized with linear optics: phase 
plates and rotators. 

Although a precise defintion of polarized light at the 
quantum level may be controversial ||, there is a wide 
consensus in viewing unpolarized states as the only 
ones that remain invariant under any polarization trans- 
formation. It turns out that the density operator of these 
states can be written as 
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\N,k)(N,k\ 



(3) 



where pn is the two-mode photon-number distribution. 
Therefore, a appears as diagonal in every subspace, with 
coefficients that guarantee the unit trace condition. 

According to our previous discussion, it seems nat- 
ural to quantify the degree of polarization of a state 
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represented by the density operator g as P(g) oc 
inf^g^ D(g, a), where U denotes the convex set of unpo- 
larized states of the form © and D(g,a) is any measure 
of distance between the two density matrices g and a. 
The degree F(g) must be normalized to the unity and sat- 
isfy some requirements motivated by both physical and 
mathematical concerns Q • There are numerous nontriv- 
ial choices for D(g,a): we shall consider here the Hilbert- 
Schmidt metric Dn S (g,(j) = \\g — a\\\ = Tr[(g — a) 2 ] 
because is the simplest one for computational purposes 
and allows us to get an analytic form of the degree of 
polarization 



Tr(g 2 



\- Pn 



^ N+l' 

N=0 



(4) 



which is determined by the purity < Tr(g 2 ) < 1 and 
the photon- number distribution pn- 

It is clear from Eq. (@J that for the states living in 
the manifold with exactly N photons, the optimum is 
reached for pure states [for which Tr(g 2 ) = 1]. In fact, 
all such pure states have the same degree of polarization: 



N 



N = 



N + l 



1 



1 

N 



(5) 



where the last expression, showing the typical scaling 
TV -1 , holds when N 1. In particular, the important 
SU(2) coherent states on the Poincare sphere [lflf 



\N,0,<t>) = ^C k (6,<f>) \N,k), 
with coefficients 



N 
fc=0 



N 



1/2 



sin- 



N-k 



cos - | - 



(6) 



(7) 



are such iV-photon states. 

However, this scaling law TV -1 can be easily sur- 
passed. Perhaps the simplest example is when both 
modes are in (quadrature) coherent states, we denote 
by | Off, ay). By reparametrizing the amplitudes as 
a H = e-'^/ 2 \/iVsin((V2) and a v = e^^yfN cos(0/2), 
we can express them as 



as, a v ) 



£< 

N=0 



- N N / 2 

=-^ 2 ^— \N,e,<t>), 



(8) 



where N = \au\ 2 + \ctv\ 2 is the average number of pho- 
tons. Since pn is a Poissonian with mean value N, one 
can perform the sum in with the result 



coh 



1 



h(2N) _ 2fl 
N 



1 



1 



2y/^N 3 / 2 



(9) 



where Ii(z) is the modified Bessel function and the last 
equality is valid for N ^> 1. 



In consequence, we are led to find optimum states for a 
fixed average number of photons N. Obtaining the whole 
optimum distribution p^ in (0J is exceedingly difficult, 
since it involves optimizing over an infinite number of 
variables. Our strategy to attack this problem is to trun- 
cate the Hilbert space and consider only photon numbers 
up to some value D, where we take the limit D — * oo at 
the end. In this truncated space, we need to find the 
states that maximize (@J with the constraints 
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o, 5>" = 1 > J2 n pk = r - ( 10 ) 



N=0 



N=0 



It is clear that the optimum must be again pure states. 
If we introduce the notations p T — (po,Pi, ■ ■ ■ ,Pd) an d 
H = 2 diag[l, 1/2, .. . ,1/(D + 1)], the task can be thus 
recast as 



1 T 

minimize — p Hp 



subject to Ap = b, 
P>0, 



(11) 



where 



b = 



A = 



1 1 1 

1 2 



1 

D 



(12) 



We deal then with a quadratic programming problem 
that, in addition, is convex, because H is positive def- 
inite |ll| . The optimum point exists and it is unique: in 
fact, there are numerous algorithms that compute this 
optimum in a quite efficient manner. Alternatively, we 
may try to determine it analytically by incorporating the 
constraints by the method of Lagrange multipliers. The 
functional to be minimized is 



£(p,\) = -p T Hp-\ T (A P -b). 



(13) 



The first-order optimality conditions V£(p, A) = to- 
gether with the initial equality constraint, give the sys- 
tem of linear equations 



H 

A 



-A T 




whose formal solution is 



A = (AH^A^-H , p = H- l A r \ 



-i AT, 



(14) 



(15) 



Before working out the analytical form of (T5J, in Fig. 1 
we have plotted the numerical solution of the quadratic 
program (|11|) for some values in < N < 1, using the 
MINQ code implemented in Matlab. The number of 
nonzero components of pn is [2A^ + 1], where the brack- 
ets denote the integer part. The distribution presents 
a clear skewness and one can check that it can be well 
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FIG. 1: Optimum distribution pjv> obtained by solving nu- 
merically the quadratic program 1111 . plotted as a function 
of the average number of photons N and N. We have taken N 
running from 0.2 to 1 and the dimension of the space D = 4. 



fitted to a Poisson distribution, which in physical terms 
means that, in this range, a quadrature coherent state 
\och, O-v) can be considered as optimum. To better assess 
this behavior, we have calculated the associated Mandel 
Q parameter 



Q 



(my 

(N) 



(16) 



where ((AN) 2 ) is the variance, which is a standard mea- 
sure of the deviation from the Poisson statistics. In Fig. 2 
we have represented Q in terms if N. As we can see, Q 
increases linearly with N and is zero only near N ~ 3. 

In Fig. 3 we have plotted the optimum distribution 
Pn for different integer values of N running from 1 to 
9. The truncation value has been chosen to be 25 in 
all the cases, although it is sufficient to ensure, for each 
value of N, that pm — for N > D. Three distinctive 
features can be immediately discerned: the solutions are 
symmetric around N, they are parabolic, and extend in 
a range from to 2N. The two first facts are in agree- 
ment with the symmetry properties of the original prob- 
lem (|llf) . The third one means a variance that scales as 
iV 2 , at difference of what happens for standard coherent 
optical processes presenting a variance linear with N (as 
for, e.g., in Poissonian or Gaussian statistics). In other 
words, the optimum states are extremely noisy and fluc- 
tuating. When N is not integer (or semi-integer), one 
can appreciate a small asymmetry that is less and less 
noticeable as N increases. 

We conclude that we can take the dimension D to be 
2N. Given the very simple form of H and A, we can 




FIG. 2: Plot of the Mandel Q parameter for the optimum 
distribution pjv obtained numerically from 1 1 I ti in terms of 
the average number of photons N of the state. 



express the final solution l|15|l in a closed analytic form: 



Pn = 3 



(N + I) 2 - (N - N) 2 3 
(2N + 1)(N + 1)(2N + 3) ~ 2iV 



N 
N 



1N*_ 
2N 2 



(17) 

which is properly normalized and shows all the afore- 
mentioned characteristics, with a maximum value of 
p N ~ 3/ (42V). If we use x = N/(2N), which can be 
considered as a quasicontinuous variable < x < 1, we 
can convert (fTTjl in p(x) — (3/-ZV) x(l — x), which is the 
Beta distribution of parameters (2,2) [13q . For the solu- 




FIG. 3: Optimum distribution pjv plotted as a function of the 
average number of photons N and N. We have taken N to be 
a integer running from 1 to 9 and the dimension of the space 
D = 25. 
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tion l|17|) , the corresponding degree of polarization is 



opt 



= 1 



(27V + l)(2N + 3) 



1 



47V 2 



(18) 



This provides a full characterization of the optimum 
states we were looking for. However, their physical im- 
plementation stands as a serious problem. The crucial 
issue for the scaling in i|18[l is the fact that distribution 
variance is proportional to N 2 . It turns out that, for the 
discrete uniparametric distributions usually encountered 
in physics, this is distinctive of the thermal (or geomet- 
ric) distribution 



1 



Pn 



N + l \N+1 



N 



N 



(19) 



But this is the photon statistics associated with the states 



10 ^7 J2 tanh " £ ^ H ® \ n ) v ' 

^ n=0 



(20) 



which are the twin beams generated in an optical para- 
metric amplifier with a vacuum-state input, with N — 
2sinh 2 £. The distribution (|19fl presents a skewness ab- 
sent in the exact solution i|17|) , but a calculation of the 
state degree of polarization gives 



1 



ln(iV/2) 
N 2 



(21) 



Apparently this is different from (|18|l . but as soon as 
N ^S> 1 they both approach unity in essentially the 
same way, which means that the (maximally entangled) 
squeezed vacuum ()20[1 is very close to optimum when 
JY > 1. In fact, one can calculate the degree of po- 
larization for (|20[1 using other metrics than the Hilbert- 
Schmidt. For example, if one employs the Bures (or fi- 
delity) distance, a simple exercise shows that p^ ures = 
1 — 1/ iV 2 , confirming again the fundamental scaling N~~ 2 . 

Before ending, two important remarks seem in order. 
First, we observe that in classical optics fully polarized 
fields have a perfectly defined relative phase between H- 
and ^-polarized modes |14| . Such a relation does not nec- 
essarily hold in the quantum domain: while the quadra- 
ture coherent states (JBJl have a sharp relative phase, the 
twin-photon beams H20J) have an almost random relative 
phase. Second, the maximally polarized states we have 
found have a highly nonclassical behavior, even in the 
limit N ^> 1, which makes the classical limit of these 
polarized states a touchy business. 
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